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ABSTRACT 


Recent SOHO I’VCS While Light Channel (WLC) observations of the 
south polar coronal hole plumes and interplume regions produce signatures 
of quasi-periodic variations in the polarized brightness (pB) at a heliocentric 
distance of 1.9 solar radii (/£q). The Fourier power spectrum of the pB time 
series shows significant peaks at about 1. 6-2.5 mHz and additional smaller 
peaks at longer and shorter time scales. Wavelet analysis of the pB time series 
shows that the coherence time of the fluctuations is about 30 minutes. The 
new observations strongly suggest that the fluctuations are compressional wave 
packets propagating in the coronal hole high above the limb. The presence of 
compressional waves may have important implications that help to explain the 
heating of coronal holes and the fast solar wind acceleration. 

Subject headings: solar wind — Sun: corona — Sun: magnetic 
fields — waves M 1 1 D 



1. Introduction 


The detection of density fluctuations in coronal holes is a difficult task due to the 
low coronal hole density (~ 5 x 10' cm”'* at 1.16 /?q, Fisher & Guhathakurta 1995) and 
the corresponding low light emission of the coronal hole plasma. Recent high-cadence 
observations by Ofrrian et al. (1997) of coronal holes using the white light channel (YVLC) 
of the Ultraviolet Coronagraph Spectrometer (UVCS) (Kohl et ah 1995) were the first to 
show polarized brightness (pB) fluctuations on a time scale of ~ 10 min high above the 
limb at 1.9 /£q. Ofmau et al. (1997, 1998) interpreted these quasi-periodic fluctuations 
as density fluctuations due to cornpressional waves propagating in the coronal hole. 

The Extreme Ultraviolet. Telescope FIT (Delaboudiniere et al. 1995) on board SOHO 
reports cornpressional slow waves in coronal hole plumes below 1.2 Rq (DeForest et al. 
1997, DeForest and Gurnian 1998, Ofman, Nakariakov & DeForest 1999). In particular, 
DeForest and Gurman (1998) observed quasiperiodic fluctuations in the brightness of Fe IX 
and Fe X line emissions at 171 A on a timescale of a few to ten minutes that coherently 
propagate outward at speeds of ~75-150 km s” 1 . Ofman, Nakariakov, & DeForest (1999) 
have identified these propagating fluctuations as slow magnetosonic waves. It is likely that 
the waves detected by Ofman et al. (1997) at 1.9 /{© are of the same nature as the waves 
detected by DeForest and Gurman (1998) lower in coronal holes. 

The energy flux contained in the slow magnetosonic waves detected by DeForest and 
Gurman (1998) is not sufficient for coronal heating. However, only part of the wave flux 
may have been detected. The cornpressional waves may carry a significant fraction of the 
energy that is required to heal the coronal holes (Ofman, Nakariakov & DeForest 1999). 
The efficient dissipation of the slow waves by compressive viscosity close to the sun, and 
by Landau damping in t he collisionless plasma suggest that the waves may contribute to 
heating of the coronal holes in the first solar radii above the limb. This region in coronal 
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holes is difficult to heal by shear Alfven waves due to their large dissipation length, even 
when phase mixing is taken into account (Ofman & Davila 1995). The compressional waves 
may provide a signature of ihe nonlinear wave acceleration mechanism proposed by Ofman 
& Davila (1997a, 1997b. 1998) for the acceleration of the fast solar wind. 

Below, we report the results of the new UVCS WLC observations of coronal hole 
density fluctuations with longer duration, and higher cadence. In § 2 we describe the 
method of observations at one and at two radial positions, together with the error analysis. 
In § 3 we show the observational results. We have applied the fast Fourier transform 
(FFT) and the wavelet analysis to our data. We performed the cross-correlation analysis 
at two radial positions in order to determine the propagation speed of the fluctuations. 
The conclusions are in I. 'Hie observations suggest that the pB fluctuations are due to 
propagating compressional wave packets in coronal holes. 


2. Observational Method 
2.1. Observations at one radial position 

The observational method of the pB with the WLC is described in detail in Ofman et 
al. (1997). Here, we present a brief summary of the method. The WLC is a polariineter 
that measures the pB in the 950-600 nm band, over a 14 x 14 square arc-second area located 
at the center of the instantaneous UVCS field of view. The coronal pB can be directly 
related to the electron density integrated over the line of sight. Because the UVCS field of 
view can be rotated about the Sun center, the WLC can look at different position angles. 
The WLC measurements and the most recent calibration are described in Romoli et al 
(1997). 

The purpose of the measurements reported below is to investigate low frequency 
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(milii-IIertz) electron density lluct.uat.ions. The determination of the absolute magnitude 
of the electron density is left for fut ure studies. The instrumental polarized stray light, 
which is usually removed from t he pB measurement, constitutes a constant and negligible 
contribution to the total pH at. the heights of interest. The time resolution is constrained 
by the required signal to noise ratio of the pB intensity. The noise in the pB intensity is 
determined by photon statistics. The number of photons that reach the detector decreases 
with the heliocentric height of the pB measurement. In order to improve the time resolution, 
we assume that the pB polarization plane is tangent to the solar limb. This allows us to 
reduce the number of polarizer positions from three to two to obtain a pB measurement. 
The pB counts are obtained by taking the difference in counts between two successive 
positions of the polarizer. 

The line-of-sight averaging of t he white light data probably reduces the variation of 
the pB due to the cornpressional waves. This effect must be taken into account when 
determining the amplitude of the waves from the observations in future studies. Looking 
at plumes may reduce the dilution due to line-of-sight integration, since the waves may 
propagate in phase due to trapping (Ofrnan, Nakariakov & Deforest 1999), minimizing line 
of sight integration effects on the signal. When interplume regions are considered only the 
largest amplitude fluctuations in the line of sight (i.e., localized density enhancements) will 
contribute to the pB intensity, minimizing the line-of-sight dilution. 


2.2. Observations at two radial positions 

In order to determine the phase speed of the density fluctuations we performed time 
resolved observations with the WLC pointing at two radial distances from the Sun with a 
cadence of about 5 minutes. This method is based on the finite propagation velocity of the 
waves carried by the solar wind. By alternating the observation of the pB between 1.9 R Q 
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and 2.1 /?© at a rate shorter than the expected travel time of the disturbances between the 
two points (on the same magnetic structure) it may be possible to observe a correlation 
between the pB fluctuations at the two points. This correlation will allow measurement of 
the propagation speed of the disturbances by calculating the travel time that yields the best 
correlation between the two time series. If the solar wind speed is independently known, 
then this travel time ('an be used to determine the phase speed of the waves, which then 
will be compared to the model predictions. 

in particular, we have used the following technique to measure the phase speed of 
the waves. The UVCS was pointed above one of the poles, and the YVLC was observing 
along the radial direction at an area of 14x14 arc-seconds. The maximum time assigned 
to continuous observation was several hours (see Table 1). During the observing run, we 
obtained three pB measurements at 1.9 Rq followed by three pB measurements at 2.1 Rq. 
This sequence was repeated for t he duration of the observation. We assumed that the angle 
of polarization did not change so that two positions of the polarizer are sufficient to derive 
the pB. As discussed below, 30-60 second exposures for each polarizer position are enough 
to achieve accuracy better than 10% on the pB measurement. The cadence was determined 
by the combination of the exposure time and the instrument mirror pointing time. 


2.3. Error analysis 

In order to estimate the signal-to-noise ratio (S/N) in our observation we define S\ and 
S ’2 as the count rates corresponding to the two polarizer positions. The pB counts are given 
by 


pB (S) - S 2 )5t y 


( 1 ) 



( 


where 6t is the exposure time (see Table 1 below). The statistical error of the pB is 

ApB - ^(S, + S.i)St. (2) 

Therefore, the S/N is given by: 

pB _ (5 1 -S a )>/5t 

ApB VST+ s 2 ' 

Computing this S/N as a function of the exposure time in seconds for the typical number 
of counts at the observed height (1.9 /£©), we get — 1.79\/£L Therefore, an exposure 
time on the order of 30 s is enough to achieve S/N— 10. This S/N is high compared to the 
S/N that can be achiever! observing the pB, with similar exposure at 1.9 jRq, using other 
presently available instruments. Although, the S/N of the pB measured with the WLC is 
high, we must bear in mind that our signal is the pB fluctuations - a fraction of the full pB. 
This reduces the actual S/N to about 3. 


A possible source of instrumental error is the WLC mirror oscillations that may 
modulate the observed brightness, leading to false detection of pB fluctuations. The 
WLC exhibits fast exponential decay right after the end of the UVCS synoptic program 
(yet to be explained), with a settling time of about one hour and irregular low frequency 
fluctuations (timescale of the order of 0.5-1 hour). The second effect does not influence the 
pB fluctuation measurement since the frequencies that we measure are much higher. The 
exponential settling time has been taken into account in the few cases where it affects our 
data sets by removing those data points from our analysis. 

In addition to white light data, we have also obtained Ly-a data with similar cadence 
and exposure time. We have investigated the variations in Ly-cv intensity using similar 
methods as for the pB data. We found that the fluctuations in the Ly-a are within the 
magnitude expected from statistical fluctuations (i.e., S/N~l). This is due to the small 
number of photons reaching the detector, resulting in poor count statistics (order of 10 
counts near the peak of the Ly-a emission line). 
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3. Observational Results 

We have made several WLC observations in February 1997, February 1998, and 
February 1999 at 1.9/k in t he south coronal hole, with exposures in the range of 30-100 
seconds. The average pH count rate was approximately 300 counts per second. Table 1 
lists the parameters of each observation: date, starting time and length, position angle 
(counterclockwise from the north pole), heliocentric height, cadence of pB measurements, 
exposure time, the frequency of the highest peak in the power spectrum, and the type of 
coronal structure (plume or interplume). For the February 1997 and 1998 observations the 
target - plume or interplume was determined by measuring a spatial intensity profile in 
Ly-a at I.oRq at the beginning of the observation and extrapolating the observed structure 
radially outward to 1.9 It^ 


3.1. FFT analysis of pB at 1.9Rq 

In Figures 1-1 we show t he observations of pB and the power spectra. The top panel 
shows the temporal evolution of the pB. The statistical relative error in the pB is on average 
9.4%. The amplitude of the variations in pB is about three times the statistical error in 
the observations. We use t he FFT spectral data analysis of the time-series (with the DC 
component and the small linear trend removed) to determine the frequency content of the 
fluctuations in pB. 

The middle panel in Figures 1 1 shows the FFT power spectrum of the pB time 
sequence. The largest peak in Figure 1 appears at a frequency / — 1.6 mHz with spectral 
resolution of 0.08 mHz. In Figure 2, the largest peak is at 1.7 mHz with spectral resolution 
of 0.03 mHz. In Figure 3, the largest peak is at 2.2 mHz. In Figure 4, the largest peak is 
at 0.5 mHz with the second largest peak at 1.2 mHz and with spectral resolution of 0.03 
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mHz. The peaks in the power spect rum of the pB indicate a periodicity consistent with 
the presence of compressional wave packets propagating in the coronal hole. The frequency 
of the peak power agrees with the expected range of frequencies for compressional waves 
driven by photospheric mot ions. 

In order to investigate the statistical significance of the power spectrum we used 
the method developed by Welch (1967): the pB time sequence was divided into ~ 1 
hour segments. The segments were windowed using the Welch window, detrended, and 
Fourier transformed, and the spectra were averaged. The average spectrum of the one 
hour segments, with the appropriate statistical error, is shown in the bottom panel of 
Figures 14. The statistical errors were obtained by multiplying the average power spectrum 
by (1 ± l/y/N), where ;V is the number of one hour segments (see Table 1). For example, 
for N = 9 the standard deviation of the average spectrum is 33%. In Figure 1, it is evident 
that there is a peak in the power spectrum at 1.5 mHz, close to the frequency of the largest 
power in the power spectrum of the full observing sequence (middle panel). The spectral 
resolution of this power spectrum is about 0.3 mHz. The spectral resolution is reduced by 
the shorter duration of the data segments compared to the full data set. The close relation 
between the location of the highest peak in the full power spectrum and the averaged power 
spectrum is also evident in Figures 2-4. The length of each observing run determines the 
standard deviation of the average spectrum. 

In Table 1, we show the frequency of the highest peak of the power spectrum that 
appears in the spectrum of each full time sequence, consistent with a peak in the average 
power spectrum of the segments at each observing run. Most of the observations exhibit 
the highest peak in the 1.6- 1.8 mHz range. In the February 4, 1998 data, the highest peak 
occurred at 0.28 mHz with a smaller peak at 2.5 mHz. The peak at 2.5 mHz appears both 
in the raw power spectrum and in the averaged power spectrum, while the low frequency 
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power of the raw power spectrum appeared spread over 1 mHz in the averaged power 
spectrum. In the February 27, 1999 data, the highest peak occurs at 2.2 mHz, close to 
the range of frequencies in the February 1997 and 1998 observations. In the lower cadence 
February 28, 1999 dat a, I he highest peak is near 0.5 mHz with the second highest peak at 
1.4 mHz. The observations at two positions (1.9 and 2.1 Rq) taken on February 5, 1998 and 
February 22, 1999 listed in 'Table 1 are described in detail in § 3.3 below. 

A way to test whether t he peaks in the power spectrum are due to aLiasing is to obtain 
spectra with different cadences. 'The observations with cadence of 74.9 s (2/7/98; 2/6/98) 
and of 89.9 s (2/24/97; 2/3/98) show the largest peak at nearly the same frequency (1.6-1. 8 
mHz). The higher cadence (54 s) observations at 2/27/99 show a peak at 2.2mHz, and the 
lower cadence (123 s) observations on 2/28/99 show the second largest peak at 1.4 mHz. To 
further test for aliasing we have decreased the cadence by a factor of two by rebinning the 
data. The power spectrum of the rebinned data still shows the largest peaks at nearly the 
same frequencies, indicating t hat the peaks are not due to aliasing. 


3.2. Wavelet analysis of pB at 1.9 R§ 

Since the wavelet transform was introduced in the early eighties (Grossmann & 
Morlet 1984), it has been widely used in signal and image processing. Like the FFT, the 
discrete wavelet transform (DWT) is a linear operation that defines a forward and inverse 
relationship between the time-domain and the frequency-domain (called the “wavelet 
domain”) with the appropriate basis functions. For the FFT, the basis functions are sine 
and cosine. For the DWT, the basis functions are more varied and complicated, and 
are usually called “wavelets.” Also, like the FFT, the DWT is orthogonal, making many 
operations computationally efficient. The advantage of the wavelet transform over the 
FFT is that the former allows the use of time-localized basis functions. This makes the 
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DWT better applicable for analysis of intermittent signals. In our case, the fluctuations 
in the pB appear as short duration (compared to the duration of the observation) pulses 
or wave-packets. Therefore, we use the Moriet wavelet (Grossmann & Morlet 1984) as our 
basis function, which is a sine wave with a Gaussian envelope 

v{t) - Tr-I^e^e-' 2 / 2 , (4; 


where the parameter % f>. 


For a given discrete time series such as our pB data, the DWT is 

\k'-k)6t] 


\v k (*) x r k ,r 

k f - o 


( 5 ) 


where s is the time scale (equivalent to wavelength in the FFT), k is the time index of the 
data points, K is the number of data points, and the asterisk denotes complex conjugate. 
The two-dimensional phase plots in Figures 5-8 below are constructed by plotting the 
wavelet amplitude for each inverse time scale 1/s (i.e. “frequency”) as a function of time 
(i.e. kSt) for the pB data shown in Figures 1 -4. 


In Figures 5-8, we show the wavelet phase diagram in which the horizontal axis is the 
time of observation, the vertical axis is the frequency in millihertz, and the colors correspond 
to the magnitude of the power at each frequency in arbitrary units. Examining the wavelet 
transforms of the pB time series it becomes evident that the quasi-periodic density 
fluctuations appear intermittently during the observations with a coherence time of about 
30 minutes. The line plot on the left panel in Figures 5-8 represents a the time-averaged 
wavelet spectrum. The good agreement between the average wavelet spectrum and the FFT 
power spectrum at the bottom panel in Figures 1-4 is encouraging and provides additional 
credibility to the spectral structure of the pB signal. 
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3.3. Cross-correlation analysis at 1.9 and 2.1 /?© 

In order to establish the wave origin of these fluctuations with higher confidence we 
have attempted to determine their group speed by observing the pB at two separate radially 
aligned points at 1.9/C and at 2.1 //©. The observations were taken on February 5, 1998 at 
14:46UT with a cadence of 89.5 s and an exposure of 60s for —9 hours and on February 22, 
1999 at 16:56LT (Fig. 9) wit h a cadence of 74 s and an exposure of 45 s for ~8 hours with 
about 4 hours of a continuous sequence (see Table 1). The WLC pointing was alternated 
between the two locations once every' 3 pB measurements. The white light brightness at 
1.9 /?© and 2.1/?© are shown at the top panel of Figure 9. As expected, the white light 
brightness of the corona is lower at 2.1/?© than at 1.9/?©. The pB at 1.9/?© (triangles) and 
at 2.1/?© (stars) are shown in t he middle panel. 

Since we expect the highest correlation between the two quasi-periodic signals when 
the fluctuations at 1.9/?© reach 2.1/?©, we look for peaks of the cross-correlation function. 
The peaks in the correlation will occur when a wave packet detected at 1.9/?© will reach 
2.1/?©. For a truly periodic wave t he peaks in the cross-correlation function occur when the 
two pB signals are in phase (i.e., shifted by a full wavelength). Negative correlation will 
occur when the waves are out of phase. A high correlation at a negative time lag may also 
occur for periodic propagating waves. 

The cross-correlation of the pB at 1.9/?© and 2.1/?© is shown in the bottom panel of 
Figure 9. In order to determine how significant the cross-correlation is - we show the level of 
cross-correlation (dashed lines) that would have been produced by two white-noise signals. 
Approximately 95% of the noise cross-correlation should lie between these bounds given by 
±1.96 /\A*i where n is the number of data pairs (Brockwell & Davis 1987). No significant 
correlation was found in the February 5, 1998 data set. However, the February 22, 1999 data 
(Figure 9) show that the cross-correlation at time lag of 656 s± 165 s is significantly above 
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the noise level and may have been caused by a propagating wave packed between the two 
points. The t ime lag error is half t he time between adjacent time lags. By dividing the 
distance between the two points (0 .211$) by the time lag of the peak correlation we obtain 
the propagation speed of the fluctuation, ~ 210 ± 54 km s~ ] , assuming radial propagation. 

The deduced propagation speed corresponds to the group speed of the compressional 
wave plus the background solar wind speed at about 2/? 0 . If we take 150±60 km s -1 for the 
solar wind speed at 2 It- based on UVCS Doppler dimming observations in coronal holes 
(Cranmer et al. 1999) we get ~ 60 ±81 km s _1 for the wave speed. The large uncertainty in 
the wave speed does not allow the determination of whether we have detected a propagating 
wave. 

The uncertainty may be reduced if we use the solar wind speed based on mass flux 
conservation at 2/? 0 , which is 90 ±20 km s _1 (Cranmer et al. 1999). In this case, we get the 
wave speed of 120±58 km s _1 in agreement with the expected sound speed in 10 6 K plasma. 
In the linear approximation (i.e., for small wave amplitudes) the waves can be viewed as 
perturbations propagating on lop of a background time- averaged, solar wind flow, where 
the time averaging is at the wave period. The mass flux conservation represents the solution 
of the equation of continuity for the background flow, on times scales larger then the wave 
periods (i.e., with the waves averaged out). We use the approximation that on long time 
scales the mass flux is conserved in order to estimate the background solar wind flow. 

In order to improve the determination of the wave speed we need to reduce the possible 
noise level by increasing t he number of data pairs. In addition, we found from the wavelet 
analysis of the fluctuations that the coherence length is of the order of 30 min. Therefore, 
we need higher cadence observat ions at the two radial distances to determine the wave 
speed with higher statistical confidence. 
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4. Conclusions 

Observations with the l VCS \\ LC channel of pB fluctuations indicate that the density 
in the south polar coronal hole at 1 .9 /{q fluctuates on a time scale of about 1.2-2. 5 mHz, 
with possibly lower and higher frequencies in the millihertz range. This is the first detailed 
high-cadence observation of such fluctuations at a considerable distance above the solar 
limb. The temporal evolution and the spectrum suggest that these fluctuations may be 
generated by compressional waves propagating from the Sun. The fact that the density 
fluctuations spectrum is peaked in a narrow frequency band at the expected range of 
frequencies strongly suggest the likelihood of wave related phenomena. Our observations 
are consistent with the recently reported (DeForest and Gurman 1998) detection of 
compressional waves by HIT in polar plumes closer to the limb (< 1.2/?q). These waves 
were identified as slow rnagnetosonic waves by Ofman, Nakariakov, & DeForest (1999). 

The fluctuations may be generated by the nonlinear compressive waves driven locally 
by nonlinear Alfven waves (Ofman & Davila 1997a, 1997b, 1998). The nonlinear Alfven 
waves provide an additional source of energy for the acceleration of the fast solar wind. The 
spatial and temporal variations of the Alfven wave amplitude in coronal holes generates 
variations in the wave pressure that drives the compressional nonlinear waves. However, 
further studies are needed to relate quantitatively the results of present observations to 
these waves. For example, the difference between the nonlinear compressional waves speed 
and the linear slow rnagnetosonic waves speed is to small to distinguish with the present 
accuracy of measurement. 

These waves could also be driven by density fluctuations at the base of the coronal hole, 
driven by the solar p-rnode oscillations that propagate into the corona in the form of fast or 
slow rnagnetosonic waves. Since the peak power of the photospheric motions occurs at 3.3 
mHz this suggests that the p-mode energy may contribute to the observed waves at about 
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half the above frequency due to nonlinear period doubling that occurs in the transition 
region. However, it is not obvious how the p-mode fluctuation could transfer their energy 
to the slow and fast mode in t he corona, and we do not attempt to solve this problem in 
the present paper. Since the energy flux in slow waves is only a small fraction of the energy 
contained in p-mode oscillations - even small conversion efficiency could potentially account 
for the observed waves. It is possible that other mechanism (such as reconnection in the 
magnetic network at the base of the corona) may generate the observed waves. 

The temporal evolution of the frequency spectrum of the pB as seen in the wavelet 
analysis indicates that the waves are produced in short bursts of about 30 minutes 
throughout the observation. The waves appear intermittently both in plume and in 
interplume regions. However, due to line-of sight effects, contributions from both regions 
are most probably contained in any single observation. 

The measurement of the propagation speed of the fluctuations produced velocities in 
the range of 150-260 km s” 1 at 2R§. This range of speed is consistent with the phase speed 
of slow magnetosonic waves or nonlinear solitary-like waves carried by the solar wind in the 
low-/? coronal hole plasma and is an encouraging result. The nonlinear solitary-like wave 
speed is about 20%-30% faster than the sound speed for typical coronal hole parameters 
(Ofrnan & Davila 1997a). Linear slow magnetosonic waves propagate at a speed close to 
the sound speed in the low-/? coronal hole plasma. Measurements of the solar wind speed 
and of the ambient temperature (in order to determine the sound speed) close to the sun 
with an accuracy of 10% or better are needed in order to distinguish between the above 
wave modes. 

The likely detection of compressional waves in the south polar coronal hole far above 
the limb may be the source of the additional momentum and heat required to accelerate 
the fast solar wind in coronal holes. The slow waves are damped by compressive viscosity 
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dose to the Sun and by l.andau damping in collisionless plasma. Preliminary calculation 
of compressive dissipat ion of slow magnetosonic waves shows that the low-frequency slow 
waves should still be apparent at 1.9 Rq for the theoretically predicted solar value of the 
compressive viscosity. It is evident that additional observations and models are needed in 
order to establish the role of these waves in the heating and the acceleration of the solar 
wind. 

This work is supported by the National Aeronautics and Space Administration under 
grant NAG5-3192 to the Smithsonian Astrophysical Observatory, by Agenzia Spaziale 
Italiana, by Swiss funding t hrough ESA's PRODEX programs, national funds and the Swiss 
Federal Institute of Technology Zurich, and by NASA SOHO Guest Investigator Program 
(NASW-98024). LO would like to thank Dr. J.M. Davila for helpful discussions. We would 
also like to thank Dr. Frederic Baudin for providing us with his wavelet transform routines, 
and Ms. Xeelima Sehgal for proofreading the manuscript. 



- 17 


REFERENCES 

Brockwell, P.J., & Davis, R.A. 1987, Time Series: Theory and Methods , Springer- Verlag, 
New York 

Cranmer, S.R., el al. 1999, ApJ, oil, 181 
DeForesl C.E., el al. 1997, Sol. Phys., 175, 393. 

DeForesl, C.E., and .J.B. Currnan 1998, ApJ, 501, L217. 

Delaboudiniere J.-P., el al. 1995, Sol. Phys., 162, 291. 

Fisher, R.R. & Guhathakurta, M. 1995, ApJ, 447, L139. 

Grossmann, A., & Morlet, J. 1984, SIAM J. Math. Anal., 15, 723 
Holder, T.E., & Leer, E. 1980, J. Geophys. Res., 85, 4665 
Kohl, J.L. et al. 1995, Solar Physics, 162, 313 
Ofman, L., & Davila, J.M. 1995, J. Geophys. Res., 100, 23413. 

Ofman, L., Davila, J.M. 1997a, ApJ, 476, 357 

Ofman, L., Davila, J.M. 1997b, ApJ, 476, L51 

Ofman, L., Davila, J.M. 1998, JGR, 103, 23677 

Ofman, L., Nakariakov, V., & Deforest, C.E.: 1999, ApJ, 514, 441 

Ofman, L., Romoli, M., Polello, G., Nod G., Kohl, J.L., 1997, ApJ, 491, LI 1 1 

Ofman, L., Romoli, M., Polello, G., Nod G., Kohl, J.L., 1998, ApJ, 507, L189 


Phillips, J.L. et al. 1995, Geophys. Res. Lett., 22, 3301 



- 18 - 


Rornoli, M., Benna, Craniner, S., Eineschi, S., Gardner, L.D., Strachan, L., Kohl, J.L., 
Nod, G., 1997, 5th SOHO Workshop, ESA SP-404, pp. 633-636 

Welch, P. D. 1967, IEEE transactions, 15, 2 

Wagner, W.J. 1976, Ap.J, 206, 583 


This manuscript was prepared with the AAS macros v4.0. 



Power Power 


19 



Time [s] 




0 0001 0.002 0.003 0.004 0.005 0.006 

f [Hz] 


Fig. 1. — Time variation of the polarized brightness observed with the UVCS YVLC in the 
south polar coronal hole at 1.9 Rq on February 6, 1998 with 45 s exposure (top panel). The 
polarized brightness units are given relative to the Sun center brightness integrated over the 
WLC wavelength bandpass (450-600 nm). The power spectrum of the pB (middle panel) and 
averaged power spectrum of nine ~ l hour segments showing the istandard error interval 
(lower panel). 
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Fig. 2. — Same as Figure 1 for observations taken on February 3, 1998 with 60 s exposure 
(top panel). Initial settling time of the mirror was removed from the analysis. The power 
spectrum of the pB (middle panel) and averaged power spectrum of nine ~ 1 hour segments 
showing the istandard error interval (lower panel). 
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Fig. 3. — Same as Figure 1 for observations taken on February 27, 1999 with 30 s exposure 
(top panel). The power spectrum of the pB (middle panel) and averaged power spectrum of 
five ~ 1 hour segments showing the istandard error interval (lower panel). 
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Fig. 4. — Same as Figure 1 for observations taken on February 28, 1999 with 100 s exposure 
(top panel). The power spectrum of the pB (middle panel) and averaged power spectrum of 
six ~ 1 hour segments showing the istandard error interval (lower panel). 









10 - 



N 

X 

JS 

O 

c 

0 ) 

0 

cr 

0) 

u 



1.2xl0 6 
1 .Ox 1 0 6 
8. Ox 10 S 
6-OxlO 5 
4.0xl0 5 
2.0xl0 S 



0 


57 115 172 230 287 


Power 


Time (min) 


Fig. 7. — Wavelet analysis of the pB data taken on February 2, 1999. Note that the strongest 
fluctuations appear in short bursts. The time-averaged spectrum is qualitatively similar to 
the averaged FFT power spectrum in Figure 3. 
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Fig. 9. — The WLC observations of white light brightness at 1.9/?© and 2.1 R Q (top panel). 
The pB at 1.9/?© (triangles) and at 2.1/i© (stars) are shown in the middle panel. The cross- 
correlation of the pB at 1.9/i© and 2.1/f© (bottom panel). The cross-correlation expected 
from white noise is marked by the dashed lines. 
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Table 1: The parameters of the high cadence UVCS WLC observations of pB fluctuations 


Date 

Start 

(UT) 

Duration 

(h) 

Position 

Angle 

Distance 

(*©) 

Cadence 

00 

Exposure 

(S) 

Peak power / Location 

(mHz) 

24 Feb 1997 

21:27 

3.7 

184.1° 

1.9 

89.9 

60 

1.8 

Plume 

3 Feb 1998 

14:19 

9.9° 

180.0° 

1.9 

89.9 

60 

1.7 

Plume 

4 Feb 1998 

14:30 

10.0 

180.0° 

1.9 

89.9 

60 

2.5 

Interplume 

5 Feb 1998 

14:47 

9.1 

180.0° 

1.9 k 2.1 

89.5 6 

60 

N/A 

Interplume (?) 

6 Feb 1998 

15:55 

8.5 

173.3° 

1.9 

74.9 

45 

1.6 

Interplume 

7 Feb 1998 

15:19 

9.1 

175.7° 

1.9 

74.9 

45 

1.7 

Plume 

22 Feb 1999 

16:56 

4.0 c 

179.6° 

1.9 k 2.1 

74.7* 

45 

N/A 

N/A 

27 Feb 1999 

16:08 

4.8 

180.0° 

1.9 

54.2 

30 

2.2 

N/A 

28 Feb 1999 

11:13 

6.6 C 

180.0° 

1.9 

122.8 

100 

0.5 

N/A 


a The pB were interpolated in the 0.40 h data gap after 4.5 h of observations. 
b It takes additional ~32 s to move the mirror between the two positions. 
c Uninterrupted observing time. 




